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c/3 ' Abstract. We study a general Majorana junction, where N helical nanowires are 

connected to a common s-wave superconductor proximity-inducing Majorana bound 
states in the wires. The normal part of each wire {j — 1, . . . , N) acts as connected 
^ ■ lead, where electrons can tunnel into the respective Majorana state 7a,j. The Majorana 

states at the other end, "/bj, are coupled to each other by an arbitrary tunnel matrix. 
We examine the conditions for even-odd parity effects in the tunnel conductance for 
' O . various junction topologies. 
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QQ ' 1. Introduction 

^ ' 

Since tlieir discovery a few years ago, many fascinating phenomena have been uncovered 
in topological insulators [1] and topological superconductors [2]. A particularly simple 
yet nontrivial example is given by the Majorana fermion bound states that must be 
present near the ends of one-dimensional (ID) topological superconductor wires; for 
^ ! reviews, see Refs. |3llll|5]. Since Majorana fermions exhibit non-Abelian statistics, they 
. are under consideration as platforms for topological quantum information processing. 
An experimental realization in terms of few-channel InSb or InAs nanowires was recently 
proposed P [7] , where the conspiracy of pronounced spin-orbit coupling, strong Zeeman 
magnetic field, and proximity-induced pairing correlations inherited from an s-wave 
superconducting substrate, will induce the Majorana end states. 

Electron tunneling from a normal metal into such a ID topological superconductor 
is predicted to show a zero-bias anomaly conductance peak due to resonant Andreev 
reflection [8l|9lll0llIIlll2lll3l|ll[T5l[T6]. In essence, the Majorana state delocalizes into 
the normal metal and forms a resonance pinned to the Fermi energy. This mechanism 
establishes a perfect Andreev reflection channel, with conductance peak value G = 2e^ /h 
at temperature T = 0. Resonant Andreev reflection is in fact the most important 
coupling mechanism in such a metal- Majorana contact as long as the proximity-induced 
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gap exceeds the thermal energy, ksT, and the apphed bias voltage scale, eV. All other 
local perturbations at the contact, e.g., quasi-particle tunneling or potential scattering, 
were shown to be irrelevant compared to tunneling via the resonant Majorana state 
^17j . Experimental reports of conductance peaks with such properties in InAs or InSb 
nanowires (with proximity-induced pairing, strong spin-orbit coupling and appropriate 
Zeeman fields) have indeed appeared recently [18], [191 1201 [21] , possibly providing first 
experimental signatures for Majoranas. Other Majorana platforms have also been 
discussed in the literature; for reviews, see Refs. [31 [U [5]. 

In order to usefully employ the nonlocal information stored in these Majorana 
states, one needs to have junctions where at least three Majoranas are coupled [221123] . 
We note in passing that regular 2D Majorana networks, e.g., with many superconducting 
grains connected by nanowires, can realize exotic spin models such as Kitaev's toric 
code [211 ES]. In Majorana junctions, the coupling between Majoranas can either 
be due to direct tunnel matrix elements, which requires that nanowires are in close 
proximity to each other, or mediated through the capacitive Coulomb charging energy 
[26l[271[28l[29l[30l[3ll[33[33llMll35]. In this paper, we study the Majorana junction 
shown schematically in Fig. [H where helical nanowires — where in the simplest 
case, right-movers have (say) spin up and left-movers spin down [1] — meet on a 
superconducting electrode. Assuming that the superconductor is grounded, the coupling 
between Majorana states is modelled phenomenologically by tunnel matrix elements. 
Note that electron-electron interactions in the normal wire parts may also imply helical 
Luttinger liquid behavior [T]. Such effects can be included along the lines of Refs. [17] [35]. 
but here we focus on the simplest noninteracting case. 

The setup in Fig. [T] provokes the question whether one can observe even-odd parity 
effects in the tunneling conductance. For tunneling into a chain of Majorana states 
[T2l [36] , there indeed is a strong parity effect: The T = conductance exhibits the 
unitary zero-bias peak value 6*1(0) = 2e^ /h for tunneling into a chain with odd number 
n of coupled Majoranas, while 6*1(0) = for even n. Moreover, n — 1 zeros and n peaks 
are present in the bias-dependent T = tunneling conductance GiiV) [12]. The strong 
parity effect is explained by noting that n coupled Majoranas, for even n, correspond 
to n/2 finite-energy fermion states. The absence of a zero-energy state then implies a 
quenching of resonant Andreev refiection, and hence Gi(0) = 0. For odd n, however, 
one zero-energy Majorana state remains and allows for perfect Andreev refiection with 
Gi(0) = 2e^ /h. For suitable choices of the tunnel couplings in Eq. ([2]) below, we will 
recover these results for the chain topology in Sec. 13. 1[ At the same time, we will extend 
them by considering other topologies of the Majorana junction, see Sees. I3.2l and[4l This 
allows us to study even-odd parity effects in a systematic manner within the general 
setup of Fig. [H 

The structure of the remainder of this paper is as follows. In Sec. [21 we describe 
the low-energy model of our Majorana junction and give the formally exact solution for 
the steady-state currents fiowing through the nanowires. This solution is elaborated for 
the "chain" and "loop" topologies of the junction in Sec. [3l based on analytical results 
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Figure 1. Schematic Majorana junction of helical nanowires (black lines) connected 
to a superconducting electrode (blue) and N metal electrodes (green); here = 3. 
The Majorana states jaj and ^bj are indicated by yellow and red dots, respectively. 
Tunneling between different jbj is encoded in the antisymmetric matrix fl, while tj 
refers to the coupling between ^A,j and 7_B.j- Tunneling from the jth lead to ^A,j is 
governed by the hybridization Tj cx Xj. 

for N = 2 and = 3, combined with numerics for arbitrary A^. In Sec. IH we show 
that analytical progress is possible for additional junction topologies, where can be 
arbitrary but Tj = T and tj = t. In particular, we address the possibility of parity 
effects in the tunneling conductance for the "ladder", "tree," and "isotropic" junction 
topologies. Finally, we offer a general interpretation of our results in Sec. [51 Some 
details have been delegated to the Appendix. 

2. Multi-terminal Majorana Junction 

2.1. Model 

We consider a junction of helical nanowires (j = 1, . . . , A^) assembled on top of a 
conventional s-wave superconducting electrode, see Fig. [1] for a schematic illustration of 
this multi-terminal Majorana junction. We assume that the superconductor is grounded, 
and hence charging effects are irrelevant. By virtue of the proximity effect, each nanowire 
segment located on top of the superconductor turns into a ID topological superconductor 
|4]. The bulk-boundary correspondence then implies that there must be a single unpaired 
Majorana fermion, 'Jbj, at the end of the superconducting wire, plus another one, 'Jaj, 
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at the interface to the normal part of the wire. The latter part also acts as connected lead 
and allows to probe electron transport through the multi-terminal Majorana junction 
in Fig. [H The Majoranas correspond to self-conjugate operators, 7^^- = 7^ ^ (with 
a = A, B), subject to the anticommutation relations 

{7a,i;7a'j'} = ^aa'Sjj'. (l) 

It is convenient to combine them to A^-component objects, •ja = (7a,i, • • • , 7a,7v)"^- 
Similarly, the conventional lead fermion operators are contained in c = (ci, . . . ,cjv)^, 
with Cj = Ylk'^jkj where cjk annihilates a lead fermion of momentum k in the jth lead. 
For simplicity, we assume identical dispersion relation for all leads, and denote their 
respective chemical potentials as fij (with f^s = for the superconductor). 

Throughout this paper, we focus on the most interesting case where the proximity- 
induced pairing gap exceeds both the thermal energy, fc^T, and the applied bias energies, 
such that quasi-particle excitations above the gap are negligible. Retaining only the 
Majoranas as relevant fermionic degrees of freedom for the junction in Fig. [H the 
Hamiltonian reads 

H = Hi+ ^-iB^lB + ilAtlB + (c^ - c)A7A, (2) 

where Hi = J2f=i ^kC^jk^jk describes the leads. The "hat" notation in Eq. ^ refers 
to an iV X iV matrix structure in wire-number (j) space. Without loss of generality, the 
antisymmetric matrix f2, containing the couplings between different 7bj, is chosen real- 
valued. The same holds true for the tunnel couplings tj, connecting 7^ and 75^- in the 
jth wire, and for A-,-, connecting 7aj- to the lead fermion c^; we use t = diag(ti, . . . , ^Ar) 
and A = (Ai, . . . , Aat). Within the standard wide-band approximation for the leads [37] . 
the Xj only appear through the hybridization scales Tj = 27rz/oA|, where f = (Fi, . . . , F^v) 
and z/q = J2k ^(^fc) is the constant lead density of states. By putting some of the 
couplings {rj,tj,Qj^ji} to zero, Eq. ([2]) serves as a general model for M < N normal- 
metal terminals connected through a junction of n < 2N coupled Majorana fermions. 
The aim of this paper is to quantitatively understand the tunneling current into such a 
multi-terminal Majorana junction. 



2.2. Conductance 

The stationary current flowing in the jth nanowire towards the junction, Ij, follows 
from the relation [121 EH] 

/, = ^ 1 de F{e - Im^^,,^,(6), (3) 

where we have Fermi distributions in the leads, F(e) = tanh(e/2A;BT), and Qaj,a'j'i^) 
the energy representation of the retarded Majorana Green's function. The main quantity 
of interest in experiments is the corresponding differential conductance Gj = edlj/dfij, 
which yields from Eq. ([3]) the T = result 

= ¥r,[-Im^:,^.^^.(/i,)]. (4) 
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Figure 2. Different topologies of a Majorana junction, shown here for = 4. Yellow 
(red) dots refer to ilBj) Majorana states, with the respective top left state for 
j = 1. A line connecting two dots indicates that the respective coupling tj or ^Ij^k is 
finite. 

Since we have a noninteracting Hamiltonian in Eq. (jj]), the Dyson equation for can 
be solved exactly. As sketched in the Appendix, the 2N x 2N matrix for thereby 
takes the form (the 2x2 block structure refers to a = A, B space) 

where 1 is the N x N unit matrix in wire- number space. The matrix element Q^jAji^) 
obtained from Eq. ([5]) then determines the current Ij and hence the conductance Gj. 

Current conservation, Ij = 0, is not imposed here due to the presence of a 

grounded superconductor [8]. Since the retarded Green's function in Eq. ([5]) does not 
depend on the chemical potentials, the current Ij, and hence also the conductance Gj, 
depends on the respective chemical potential fij only, but not on the fj.k^j. All "nonlocal" 
conductances, edlj/d/ikj^j, therefore vanish identically regardless of the values of the 
tunnel couplings, as long as the validity requirements for Eq. ([2]) are met. The resulting 
decoupling of different terminals is a characteristic feature of noninteracting Majorana 
networks and can be traced back to the assumption of a grounded superconductor. 
For a "floating" (not grounded) mesoscopic superconductor, where current must be 
conserved and Coulomb charging effects may become important, nonlocal conductances 
are typically finite. 

From now on, without loss of generality, we consider only the tunneling conductance 
Giiy) = dIi{V)/dV for nanowire j = 1, with eV = Hi and in the most interesting T = 
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limit, see Eq. (jl]). It is worth stressing that while Gi is independent of the chemical 
potentials /ij>i, it still does depend on all hybridization parameters Tj. Although Gi{V) 
can always be expressed in terms of the retarded Green's function ([5]), explicit results 
for arbitrary parameters are generally lengthy and not illuminating. We here instead 
aim at analytical results by examining several limiting cases for the specific junction 
topologies drawn in Fig. [2j In Sec. [3l we consider N = 2 nanowires with arbitrary 
coupling matrix elements, where a compact expression for Gi{V) allows us to study 
tunneling into a "chain" of Majoranas. The case of = 3 wires is also addressed, 
where we obtain insights for the "loop" geometry in Fig. [2l We then analyze the linear 
conductance Gi(0) for loop configurations with arbitrary N > 2 Majorana states in the 
loop and study the parity effect for such a junction topology. In Sec. HI we discuss the 
remaining three topologies shown in Fig. [21 where we assume Tj = T and tj = t in order 
to simplify the analysis. In addition, all non-zero matrix entries in fl are assumed equal, 
i.e., flj<:k is either some finite constant (Q) or zero. Under these conditions, we show 
that analytical results for Gi{V) with arbitrary can be derived, which allows us to 
systematically study even-odd parity effects in Majorana junctions. Before moving on, 
however, we first discuss a useful reformulation of Eq. ([5]). 



2. 3. Retarded Green's function 

To achieve analytical progress for arbitrary A^, at least for some special parameter 
choices in Eq. ([2]), it is convenient to write the relevant matrix element, Q^^^ ^1(^)5 '^^ 
retarded Green's function in an alternative form avoiding direct matrix inversion. We 
first single out the effects of the inter-wire couplings Cl by writing a Dyson equation for 

g^-Go^Go^l ;^)^^ (6) 

where the "unperturbed" {Cl = 0) retarded Green's function Qq is diagonal in wire- 
number space. In fact, we can read off its matrix elements from Eq. ([5]), 



e 



a,a' 



Now we introduce the auxiliary quantities 

which encode the off-diagonal matrix elements of the retarded Green's function and 
fulfill the relation 

X, - ^ ' ' ^ E n,,X, = 5,,,. (9) 
e + j k=i 
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Equation follows by taking the {Bj,Al) matrix element of Eq. (jS]), and is widely 
used below to determine the Xj. Note that for the tunneling conductance (jl]) we need 
GaiaiI^)- Taking the (^41,^41) matrix element of Eq. ([6]) yields after some algebra 

Equation (JTOll can be further simplified for e = eV — > 0, where the linear conductance 
is 

Gi(0) = ^(l-Xi). (11) 
The real- valued Xi = Xi(e = 0) follows from solving Eq. (Q. 

3. Chain and loop junctions 

3.1. Conductance for chain topology 

Let us start with the case of = 2 nanowires, where Cl is fully determined by just 
one parameter, Q12 = ^l. Equation ([5]) then yields the Green's function matrix element 
determining the tunneling conductance, 

,x etl + in^~e^)ie + tT2) 



/Ai,AU ) ^2(, + ,ri)(e + tV^) - (e2 -tl + ^eT^){e^ - + ^eEa) 



(12) 



Several remarks about this result are in order, (i) For ti = 0, Eq. ( |T2l) reduces to 
Qaiai — V(^ + ^^i), which leads to the well-known resonant Andreev reflection peak 
for tunneling into a single decoupled Majorana fermion (7A1) [8l [TT]. 

Note that the parameters ^2 and T2 do not affect this result at all. (ii) Without coupling 
between the wires, f2 = 0, but with ti 7^ 0, the Majorana jb,! is also involved and the 
resonance ( !T3|) will split into two symmetric finite-voltage peaks separated by a perfect 
dip at zero bias [U [TT]. Indeed, now Eq. ( fT2l) yields 

2e2 (el/El 

T{evr,y + [{evy-tl] 

such that the linear {V — )■ 0) conductance vanishes but two perfect (Gi = 2e^/h) 
peaks are present for eV = ±ti. This reflects the parity effect in a chain of n coupled 
Majoranas, where Gi{V) has n — 1 zeroes and n unitary peaks [12], in accordance with 
Eq. ( Jl4l) for n = 2. (iii) Putting ^2 = in Eq. ( IT2l) . we effectively recover tunneling into 
the n = 3 Majorana chain composed by 7^1 —751 — 1b2^ 

2e^ Vlp?-{eVf '^ 

T {eV)^[n^ + tl~ (e\/)2]2 + Tl[n^ - {eVf 

which exhibits three unitary conductance peaks, at V = and eV = ±\/Q'^ + tf, and 
two zeroes, eV = ±Q, as expected under the parity effect, (iv) Finally, the n = 4 



) /-TM-jro^J , -i2 _ /^T/\212 I r2r02 _ /'^T/^212' 
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Majorana chain is realized for r2 = in Eq. fll2p . For simplicity, we put ti = t2 = t and 
obtain 

^ ^ 2e^ {eVT,nn' + t'-{eV)r .^g. 

' h [{evyQ^-{{evy-t^yf + {evriy[Q + t^-{evy]^' ^ ' 



Clearly, we have three zeroes, = and eV = iyHM-t^, and four unitary conductance 
peaks, eV = ±^ ± a/T^^/T+T^ with independent ± signs. When also r2 7^ 0, the parity 
effect is not ideal anymore, i.e. peak conductances are smaller than 2e^ /h, but it remains 
observable. 

The above Majorana chain configuration can be generalized to setups with N > 2 
and n < N + 1 Majorana states, see Fig. [2l Let us focus on computing the linear 
conductance G'i(O) for n = + 1. To that end, we put tj^i = and keep as non-zero 
entries in Cl only the nearest-neighbor couplings ^jj+i with j = 1, . . . , — 1. In the 
limit e — )■ 0, the Dyson equation ([9]) yields 

T1Q12 

Xi 72 — X2 = 1, Qjj+iXj^i + Qjj_iXj_i = (17) 

^1 

for j = 2, . . . , N. The boundary condition X^r+i = leads to Xj odd = {Xj even = 0) 
for even (odd) A^, and hence Xi = S^^odd- The conductance (ITT!) for tunneling into this 
chain of n = + 1 Majoranas is 

Gi(0) = ^(5„,odd, (18) 

and exhibits a strong parity effect. Moreover, we have numerically confirmed the strong 
parity effect in such a junction by computing the differential conductance. In full 
agreement with Ref. |T2], we found that Gi{V) exhibits n unitary-limit (2e^//i) peaks 
as a function of bias V. 

3.2. Conductance for loop topology 

The algebra needed to obtain Gi{y) from Eq. (j5]) for A^ > 2 Majorana wires becomes 
much more involved. For A^ = 3, the linear conductance can still be given in compact 
form for arbitrary parameters. We find: 

2e^ Fif 

'''<°' = -f^FTifSPlWsiJ ■ i'^^^^'^A- (19) 

For ti 7^ 0, the unitary limit {Gi = 2e^/h) is reached when either t2 = ^3^23 = or 
= T2^23 = 0, while at the same time ^ 7^ 0. Note that 6*1(0) = when ^2 = ^3 = 0. 
The junction then consists of 7^1 coupled to a closed "loop" consisting of the three B 
Majoranas. This is the simplest example for the "loop" topology shown in Fig. O 

The linear conductance Gi(0) for loop junction configurations with arbitrary N > 2 
and non- vanishing fljj+i can readily be analyzed using the Dyson equation ([9]). As 
before, see Sec. 13.11 assuming tj = and taking the limit e — )■ in Eq. ([9]), we obtain 

Xi — -J (^12X2 + QinXm) = 1, fijj+iXj+i + Qjj^iXj^i = (20) 




Figure 3. Tunneling conductance Gi{V) vs. eV/T (witli F = Fi) for a Majorana 
junction in the homogeneous "loop" configuration {ftjj^i — ft), shown for several N 

with ti/r = n/r = i. 



for j = 2, . . . , N , with periodic boundary conditions, X^v+i = Xi. For odd A^, one can 
verify that Eq. f l20p has always, i.e., for arbitrary Qjj^i, the solution Xi = 1. Indeed, 
for Xi = 1, the two sequences Xj odd/even, both obeying 

^. = %^^.-2 (21) 

with j 7^ l(modA^), are uniquely determined by virtue of the matching relation 
0,12X2 + QinXn = 0. As a result, the linear conductance f|TT]) always vanishes for 
odd A^. For even A^, a general solution to Eq. fl20|) is given by 

Xj odd =0, (22) 



-^j even ^ ^j—2i 



^^12X2 + ^IinXn = —tl/Ti, 

implying that the linear conductance f lTTj) reaches the unitary-limit value 2e^ /h. This 
suggests that like for a Majorana chain, the strong parity effect described by Eq. f|T8l) 
with n = A^ + 1 is also present for the loop topology. Thus, by changing the number of 
tunnel coupled nanowires A^ in the loop configuration, one may switch the conductance 
Gi (0) on and off, depending on the parity of A^. This parity-based switching mechanism 
for the tunneling conductance stems from the non-local nature of electron transport in 
Majorana junctions, and offers a way to "engineer" such parity effects. 




Figure 4. Same as in Fig. [3] but for generic inhomogeneous "loop" configuration, 
wliere only ili^2 — O.GSFi has been changed; all other ^i^i+i = Ti as in Fig. [S] 

However, the solution fl2^ fails when is even and the parameters ^jj+i satisfy 
the incompatibility condition 

+ =0. (23) 

, "7 "2fc+l,2fc+2 

In particular, for a homogeneous loop configuration with identical couplings [^jj+i = ^ 
for all j], Eq. (123|) is satisfied. We then have the trivial solution Xj = 1 for all j, 
regardless of the parity of > 2. Consequently, Gi(0) = for arbitrary (odd or 
even) number of Majoranas in such a loop. This can be understood as a result of the 
destructive interference - with complete cancellation of clockwise and anti-clockwise 
contributions to 6*1(0) - in a homogeneous Majorana loop (closed chain). 

Remarkably, for even and arbitrary Majorana loops, this interference therefore 
results in either a maximal {Gi = 2e^//i) or a completely blocked {Gi = 0) zero- 
bias conductance, depending on the incompatibility condition in Eq. f l23|) . Fig. [3] 
shows the differential conductance Gi{V) for several in the homogeneous loop 
configuration. Numerically, we find that the number P of finite-bias peaks in Gi{V) for 
this configuration is given by 

N + 1, AT odd, 

2{l + [N/4\), AT even, ^ ' 

with [x] denoting the integer part of x. Due to the lack of an unambiguous 
correspondence between P and A^, a parity effect is difficult to detect in the homogeneous 
loop. 
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Let us stress again that in the loop configuration, Gi{V) exhibits the strong parity 
effect, with + 1 unitary-hmit (2e^//i) peaks as in the chain configuration, for generic 
(inhomogeneous) parameters. This is illustrated in Fig. H] for a particular choice of the 
couplings. 



4. Junction topologies and parity effects 

In this section, we address the three remaining junction topologies in Fig.[2]for arbitrary 
A^. The tunneling conductance is obtained for a parameter set with equal hybridizations, 
Tj = r, and intra-wire couplings, tj = t. The matrix elements Qj^k are either zero or 
equal to a single constant (fi), depending on the junction type. 



4.1. Ladder configuration 

For the "ladder" shown in Fig. [21 the only non-zero fij<:k entries are nearest-neighbor 
bonds, ^j,j+i = fl for j = 1,...,N ("closed ladder") or j = 1,...,A^ — 1 ("open 
ladder"). The example shown in Fig. [prefers to an open ladder for = 4. For both 
ladder configurations, Eq. ([9]) for Xj{e) takes the form 

X,-A(X,+i-X,_0 = 5„i, M^) = (25) 

supplemented by the boundary conditions Xq = X^+i = {Xj = Xjj^n) for the open 
(closed) ladder configuration. 

Starting with the open ladder, the solution of Eq. is 

v/r+4l^F^+^ + F_^+^ 1 l±x/r+4A^ 
^U^) - — ^ — - ^±iej ^ • (ibj 

The conductance Gi{V) follows readily from Eqs. (jl]) and ffTOj) . For ^ 1, using 
> 1, we find that Xi in Eq. fl2B]) becomes A^- independent, Xi ~ —F_/A. Hence 
A^-dependent parity effects must disappear in this junction type for large A^. However, 
for moderate A^, we still observe even-odd parity effects, see Fig. [5l For odd (even) A^, 
Gi{V) exhibits a dip (peak) near zero bias. However, peaks do not reach the unitary 
limit (Gi = 2e^ /h) and dips are not associated with perfect conductance zeroes anymore. 
We call this behavior "weak parity effect" to contrast it from the ideal "strong parity 
effect" found for tunneling into a Majorana chain or loop, see Sec. [31 Notice that the 
even-odd features, clearly visible in the lower inset in Fig. [5l for the linear conductance, 
gradually disappear as A^ increases, as expected from the above argument. The upper 
right inset of Fig. [5l shows that the linear conductance Gi (0) monotonically increases 
with the effective parameter A(0) = VtV /t^. The tunnel couplings in our model enter 
Giiff) in this configuration only through A{f)). 

For the closed ladder, we find a very similar picture. Now the solution to Eq. ([25ll 
follows by Fourier expansion, 

1 1 

^^^^^ ^ N 1 - 2iA{e) sm{7im/N) ' ^"^^^ 




Figure 5. Differential conductance Gi{V) for tunneling into a Majorana junction in 
the "open ladder" configuration. Main panel: Gi vs. eV/T with t/T = 0.5,n/r = 1.2 
and several N. Upper right inset: Linear conductance Gi(0) vs. H.T/t'^ for iV = 3,4, 
and 10. Lower inset: iV-dependence of Gi(0) for same parameters. 



In the large- limit, this expression turns into an A^- independent integral and parity 
effects are quenched. Numerical evaluation [based on Eqs. (jlj), ( ITO!) and (l27j) ] shows that 
that the differential conductance Gi{V) has P peaks, where the number P coincides with 
the one found in the homogeneous loop configuration, see Eq. fl2^ . However, now the 
peak structure is less pronounced. 

To conclude, a Majorana junction with ladder topology, either of closed or open 
type, is found to exhibit a weak even-odd parity effect that disappears for N ^ oo. 

4.2. Tree configuration 

Next we study the "tree" topology in Fig. [2|, where 7^,1 is connected to all other 7B,j>i 
through Qi j = Q, but the remaining f2-couplings vanish. Now Eq. (jO]) has the solution 

The linear conductance then follows from Eq. ( ITTj) . 



Evei 



~\ — I 1 1 1 — I — I 1 1 1 — I — I 1 1 1 1 II 



_i I I I I 




Figure 6. Gi{V) for the "isotropic" case with N — A and N — 5, for t/T = 0.5 and 
ri/r = 1.2. Inset: Linear conductance Gi(0) vs. A = QT/t'^ for iV = 5. 



There is no oscillatory (even-odd) A^-dependence, and we find no parity effect in the tree 
configuration. Interestingly, as increases, 7^,1 effectively decouples and the unitary 
single- Major ana limit is approached. 



4-3. Isotropic configuration 

Finally, consider the isotropic configuration in Fig. |2l where all couplings Qj^k = ^ 
have the same magnitude. (Note that similar isotropic couplings also arise through the 
charging interaction [5^ 155] .) As detailed in the Appendix, for the isotropic Majorana 
junction with arbitrary A^, the tunneling conductance is given by 



2e2 

G,iV) = — Im 
Xi(e) 



e + iT 
1 



1 



e2 - t2 + 



=eV 



where A{e) = iQ{e + iT)/{e^ - t^ + ieV) is as in Eq. 



6(e) 



1 + A{e] 
1 - A{e] 



N-l 



(30) 

for the ladder topology, and 

(31) 



The nonlinear tunneling conductance Gi{V) is shown in Fig.|6l It exhibits a weak parity 
effect, with a zero-bias peak (dip) in Gi{V) for even (odd) A^. In particular, the dip 



Even- odd parity effects in Majorana junctions 



14 



does not lead to a perfectly vanishing linear conductance. 

To study this issue in more detail, note that for e = eV = one has A = QT/t^ 
and the linear conductance is given by Eq. (fTTi) . For A ^ N, Eq. ( 1311) implies 
^1 ~ (^Af.odd + (— 1)^~^(A^ — 1)M, and hence we get in this limit Xi ~ d^^odd/N, i.e., 



The linear conductance thus exhibits a parity effect for QT/t'^ ^ 1. However, this effect 
disappears for large A^. The approach to Eq. fl32p with growing A is shown in the inset 



We conclude that the isotropic Majorana junction is characterized by a weak parity 
effect for small-to-intermediate A^, quite similar to the ladder configuration in Sec. 14.11 

5. Discussion and conclusions 

We now offer an intuitive and general interpretation of our results for the linear tunneling 
conductance Gi = Gi(0), where lead j = 1 is connected to a set Ai of coupled Majorana 
fermions. Every Majorana fermion in the set A4 is linked to other Majoranas in this set 
through K = 1,2, . . . "bonds" , where a bond requires that the respective tunnel coupling 
is finite. In general, the set Ai contains not only the jaj but also additional Majorana 
fermions rjj oc Cj — corresponding to the normal parts (leads). The rjj are coupled only 
to 'jAj by \j, see Eq. ([2]), with Tj oc A^. We note in passing that within the bosonization 
approach, the rjj represent Klein factors for each lead [M| [35]. A Majorana mode with 
K = 1 (but not rji) is called "outer Majorana fermion" and, if present, constitutes a 
boundary of the set Ai. The tunneling conductance Gi then follows by summing over 
all propagation amplitudes for paths on Ai connecting r/i and outer Majoranas, where 
a possible parity effect follows from the following rule: Trajectories containing an odd 
(even) number of Majoranas in the set M. yield a finite (zero) contribution to Gi. If 
has no boundary (i.e., no outer Majoranas exist), Gi is determined by the sum over all 
"loop" trajectories passing through 7^1. In particular, "homogeneous loop" trajectories 
do not contribute to Gi, see our analysis in Sec. 13.21 

We have verified this rule for the case of A^ = 3 nanowires (for arbitrary parameters) 
in Sec. 13.11 as well as for all the examples in Sec. |H For instance, in the "tree" 
configuration, Eq. (|29l) shows that Gi 7^ for arbitrary A^. Since there are A^ — 1 
outer Majoranas ?7j>i, each connected to rji through a chain of four Majoranas, every 
path connecting rji and rjj^i involves an odd number (five) of Majoranas. (Note that 
rji is not included in the set M., and hence, in the Majorana number, it should not be 
counted while the outer Majorana rjj^i counts.) As a result, a finite conductance Gi 
follows from the above rule. As another example, consider the "ladder" configuration 
with A^ wires, where we also have A^ — 1 outer Majoranas (?7j>i). Here, paths connecting 
rji and rjm-, for even m = 2, 4, . . ., lead to an enhancement of Gi, which can explain the 
reported peak vs. dip structure for even vs. odd A^. 




(32) 



of Fig. E] for A^ = 5. 
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The present work has shown that rich parity effects are present in junctions of 
multiple Majorana nanowires. For chain-like arrangements, the T = linear tunneling 
conductance is either zero (for even number of Majoranas) or 2e^ /h (for odd number) 
|12j . Such a strong parity effect does not survive for general junction topologies. 
For instance, in the "ladder" and "isotropic" topologies, see Fig. [2], one has only a 
"weak parity effect", where non-ideal zero-bias dips and peaks can be observed in the 
tunneling conductance depending on the parity of the Majorana number. However, 
these effects are less pronounced and disappear with increasing number of Majoranas. 
In other configurations, e.g., for the "tree" topology, there is no parity effect at all. A 
particularly noteworthy result concerns the tunneling conductance when '-yA,i is coupled 
to a "Majorana loop." For generic parameters, a strong parity effect is present for such 
a junction topology, with the even-odd features in the conductance as for a Majorana 
chain. However, for a homogeneous loop (and in some other special cases), the linear 
conductance is found to vanish for any number of Majorana fermions in the loop, i.e., 
no parity effect can be detected. 

Finally, let us briefly address the role of finite-temperature effects and of Coulomb 
interactions. Finite T implies a thermal broadening of all peak or dip features in the 
differential conductance, plus a reduction of the peaks. Once T reaches the energy scale 
corresponding to the separation between two adjacent peaks in GiiV), peaks start to 
overlap and the parity effects described here will be smeared out. Concerning Coulomb 
interactions, it has recently been pointed out that charging effects cause interesting 
Kondo physics in similar Majorana networks [321 EH BSlj. This physics is different from 
the parity effects studied here, and it would be interesting to study their interplay with 
Kondo physics for finite charging energy. 
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Appendix 

Retarded Majorana Green's function 

Here we briefly sketch the derivation of Eq. ([5]), obtained by solving the combined 
equations of motion for Majoranas and lead fermions, 

j' 

{idt - ik)ckj = \lA,j- 
The retarded Green's function is contained in the general Keldysh Green's function 
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where Tc denotes time ordering along the standard Keldysh contour [3Zj. Exploiting 
the wide-band approximation for the leads, some algebra yields the Dyson equation for 
Q from the above equations of motion, 





The 2x2 structure refers to a = A,B space. In addition, every matrix entry 

still carries the standard 2x2 Keldysh substructure [37], ^ = I ^ ^"1' ^^^^^ 

g°- = [g^]^ and the Keldysh component is . The Keldysh structure of the self-energy 
S(e) = diag(Si, . . . , Sat) due to the lead fermions is 



1 F(e-/x,) + F(e + /i,: 
-1 



where F{e) = tanh{e/2kBT). Noting that Cl and i are diagonal matrices in Keldysh 
space, it is then straightforward to extract the 2N x 2N matrix for the retarded Green's 
function g^ quoted in Eq. ([5]). 

Isotropic junction 

Here we show how to derive the quantity Xi{e) quoted in Eq. fl30|) for the isotropic 
topology in Sec. 14.31 First, it is useful to define auxiliary variables Yj{e) = J2i=i-^k{^) 
for j = 0, . . . , N, with Yq = 0. The Dyson equation for Xj, see Eq. (|9]), then takes the 
form 

(1 + A)Y, - (1 - = Sj,i + AYn, 

with A{e) specified in Eq. f l25|) . For j = 1, this reduces to 
^ _1 + AYn 

^'-^Ta~' 

while for j > 1, we arrive at a non- homogeneous linear recursion relation for the 
quantities = Yj/Y^, 

(i + A)e,-(i-A)e,-i = A, e^ = i. 

One can easily verify that this set of equations is solved by 
1 



2 



1 + 



1-A 



where Eq. ( 13T|) follows for j = 1. Now note that Xi = Yi = ^lY^ but also 
Yi = (1 + AYi\i)/{l + A), see above. We thereby obtain Xi(e) in the quoted form 
[Eq. (I30D], where Gi{V) follows from Eqs. and (fTOl). 
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